Step-like magnetization in a spin-chain system: Ca3Co206 
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Due to a ferromagnetic in-chain coupling between Co^^ ions at trigonal sites, chains C02O6 are 
considered as large rigid spin moments. The antiferromagnetic Ising model on the triangular lattice 
is applied to describe an interchain ordering. An evolution of metastable states in a sweeping 
magnetic field is investigated by the single-flip technique. At the first approximation two steps 
in the magnetization curve and a plateau at 1/3 of the saturation magnetization are found. Four 
steps in magnetization are determined in high-order approximations in agreement with experimental 
results. 

PACS numbers: 75.25.-|-z, 75.30.Kz, 75.50.Ee 



Among other spin-chain compounds, Ca3Co206 has 
drawn recently considerable attention to their complex 
magnetic behayior ill ^, .3, 4, 5j. The most intriguing 
feature observed in Ca^Co2p6 is a step-like shape of the 
magnetization curye [^|^,|^. The number of the steps in 
the curve depends strongly on a sweep rate of the exter- 
nal magnetic field and temperature |2, S Q ■ Two steps 
become apparent in the temperature range from 12 K to 
24 K j^. The first step takes place at the zero magnetic 
field. Then the magnetic moment remains constant at 
about 1/3 of the full magnetization up to the magnetic 
field of 3.6 T where the second step occurs to the fully 
magnetized FM state. At least four equidistant steps are 
clearly visible below 10 K at a very low sweep rate. They 
are accompanied by a sizeable hysteresis. Similar phe- 
nomena were observed in other spin-chain compounds, 
e.g. CaaCoRhOe 



The structure of Ca3Co206 consists of C02O6 chains 
running along the c axis. The Ca ions are situated be- 
tween them. The chains are made up of alternating, face- 
sharing CoOg trigonal prisms and CoOg octahedra. The 
crystalline electric field splits the energy level of Co^+ 
ions into the high-spin {S = 2) and low-spin {S = 0) 
states. The Co"^"*" ions situated in the trigonal environ- 
ment (Col) are in the high-spin state and the octahedral 
Co sites (Coll) occurs in the low-spin state. In the last 
case the energy difference between the low-spin and high- 
spin states is very small and a tiny fraction of Coll sites is 
reported to be in the high-spin state. The crystalline elec- 
tric field leads also to a very strong Ising-like anisotropy 
at the Col sites. The chains form triangular lattice in the 
ab plane that is perpendicular to the chains. An in-chain 
exchange interaction between magnetic Col ions through 
the octahedra with non-magnetic Coll ions is ferromag- 
netic (FM). It causes the in-chain FM ordering of Col 
ions at about 40 K. The interchain interaction is antifer- 
romagnetic (AFM) and much weaker than the in-chain 
one. A partial AFM order of chains appears at 24 K. A 



weak feature concerned most probably with a transition 
in a new interchain order was also observed at around 12 
K. This scenario of magnetic interactions in Ca3Co206 
is consistent with results of x-ray photoemission spec- 
troscopy Q , neutron scatteringl9|, magnetization and 
specific heat measurements 0, 0, S ^ IM nuclear 
magnetic resonance and theoretical calculations of 
indirect interactions between Col sites p^ . 

The model presented in Ref. deals with an in-chain 
structure of Ca3Co206 and magnetization dynamics ex- 
plained in terms of the quantum tunnelling. In the this 
Letter, we develop a new model for a description of the 
step-like magnetization in Ca3Co206 at low tempera- 
tures, shifting the stress on the interchain magnetic or- 
der. The strong FM in-chain coupling makes it possible 
to consider a Ca206 chain as a large rigid spin formed by 
Col ions. There are only two projections of the chain spin 
onto the c axis due to the strong Ising-like anisotropy. In- 
cluding the AFM coupling between the nearest-neighbor 
chain spins we arrive to the Ising Hamiltonian on the 
triangular lattice 



H = J 



<iO> 



(1) 



where a\ — ±1 is the c-axis projection of the i-th chain 
spin, J > is the parameter of the AFM interchain cou- 
pling, B is the magnetic field, < ij > denotes the summa- 
tion over all the nearest-neighbor pairs on the triangular 
lattice. 

The strong dependence of the magnetization curve 
shape on the magnetic field sweep rate and temperature 
shows that the state of the system of the chain spins in 
the magnetic field is far from equilibrium. At the low 
sweep rate the system is rather in a metastable state 
than in the ground state. It is convenient to formulate 
necessary conditions of the metastability of the system 
in the following form 

(2) 



a^hi < 



where 
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is the effective field for the i-th chain, denotes sum- 
mation over the nearest-neighbors of the i-th chain. We 
have used the unstrict inequahty in Eq.|21) keeping in 
mind a strong degeneracy of partially ordered AFM ar- 
rangements on the triangular lattice. A transition from 
one metastable state to another occurs through exited 
states. We obtain a-fhi > at least for one chain in a 
exited state by definition. Let AEi = 2afhi > be the 
excitation energy per Col site of the i-th chain. It fol- 
lows from this that the probability of an exited chain at 
low temperature T is extremely small oc exp{—NAEi/T) 
since the number of Col sites in the chain (iV) is consid- 
ered to be large. That is why, we should investigate an 
evolution of metastable states in the slowly sweeping ex- 
ternal magnetic field assuming T = quench. 

We perform the investigation of the evolution of the 
system using the single-flip technique that was applied 
earlier to nonequilibrium dynamics of the AFM Ising 
model on the the triangular lattice fT3|. In the terms 
of the effective field, the spin-flip probability A is taken 
in the following form 



A 



if afhi < 0, 

1 ifcr?/i,>0. 



(3) 



In contrast to Ref. where the spin-flip technique 
was applied to numerical Monte Carlo simulations, we in- 
vestigate the evolution analytically. If a state under con- 
sideration is degenerate and different sequences of spin 
flips lead to different final states, one should take into 
account each possible sequence with equal probabilities. 
This assumption can be proven rigorously by Bogolubov's 
quasiaverage technique. That is a small auxiliary random 
field should be added to the Hamiltonian ^ in order to 
lift the degeneracy. After that we take an average over 
a manifold of the auxiliary fields restoring equivalence of 
different chains and, then, let the amplitude of the auxil- 
iary fields goes to zero. This approach is equivalent to the 
Monte Carlo technique in the limit of the large number 
of samples. 

We take the ground state of the triangular lattice as 
an initial state of the chain lattice at i? = 0. The ground 
state of triangular lattice at B = is strongly degener- 
ate and it is impossible to represent it in an explicit form 
[Ti'l . On the other hand, we can produce a set of approx- 
imations for the ground state. Since the entropy density 
for the ground state was calculated exactly by Wannier 
fil as 



S=- 



ln(2cosw)dw « 0.3383, 



(4) 



we are able to compare the entropy density of the ap- 
proximated state to the exact value in order to control 
the precision of our approximation. Previously, Maignan 
et al. 1^ performed a qualitative analysis of the magneti- 
zation curve starting with an initial state that consisted 
of alternating rows of spin-up and spin-down chains. The 
energy of such an arrangement equals the ground state 
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FIG. 1: The honeycomb magnetic structure. The black and 
white circles are spin-up and spin-down states, correspond- 
ingly, the gray circles can be either spin-up or spin-down 
states. 



energy but its statistical weight goes to zero in the limit 
of the infinite lattice. Therefore, arrangements of this 
type should be discarded Q . 

The first approximation to the ground state of the 
AFM Ising model on the triangular lattice is the honey- 
comb structure shown in Fig^ Two thirds of the total 
number of the chain spins are ordered in the AFM honey- 
comb structure whereas other chain spins placed in the 
centers of hexagons have arbitrary projections of chain 
spins onto c axis. The entropy density of the first ap- 
proximation equals Si = (l/3)/n2 « 0.231. An arbirary 
small external magnetic field lifts degeneracy of chain 
spins in the centers of hexagons orienting them along the 
magnetic field (trf — 1). We consider that the spin-up 
chains (the black circles) are directed parallel to the mag- 
netic field. The grey circles become black in FigQ] This 
causes a step at _B = with the height of 1/3 where 
the full magnetization is taken as unity. Then, one third 
of the spin chains remain spin-down or antiparallel to 
the field. Since they are surrounded by 6 spin-up chains 
this configuration is stable up to the critical magnetic 
field Be = 6J where a transition to the fully-polarized 
FM state takes place. The magnetization curve for the 
first approximation is shown in Fig|21 It should be no- 
ticed that this curve is in an excellent agreement with 
the experimental data at the intermediate temperatures 

mm 

We can improve the approximation used above includ- 
ing tripod configurations |l4| . Since the chains placed in 
the centers of hexagons have an arbitrary projections of 
the spin onto c axis, it occurs that three of them neigh- 
boring with the same chain of the honeycomb sublattice 
are in the same state. Then the configuration shown in 
FigO can appear. The chain spin in the center of the 
tripod belongs to the honeycomb sublattice but, nev- 
ertheless, it can have an arbitrary c projection of the 
chain spin. These states increase the entropy density 
up to S" = (5/12)Zn2 « 0.289 ^ and involve various 
types of configurations: isolated tripods and connected 
tripods. It is convenient to consider them separately. Af- 
ter straightforward calculations we obtain the probability 
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FIG. 2: The magnetic moment as a function of the dimen- 
sionless magnetic field BjJ. The dash, dash-dot, and sohd 
hnes are the results of the first, second, and fourth approxi- 
mations, correspondingly. The four nearest-neighbor configu- 
rations producing the critical spin-fiip fields are shown. Note 
the break in the vertical axis. 



of an isolated tripod 

Pa = (1/12)(1 - 1/2^ - 1/2'' - 1/2^)3 « 0.0397. (5) 

We include this type of configurations in the second ap- 
proximation for the initial state. The entropy density 
of the second approximation is 5*2 ~ 0.259. Applying 
the external magnetic field we again obtain a step at the 
zero magnetic field but the chain spins in the centers 
of hexagons sharing joint corners with the tripod cen- 
ter remain in the spin-down state because they are sur- 
rounded by 4 spin-up and 2 spin-down nearest-neighbors. 
The height of the step (AM(i?/J)) can be expressed 
through probability of the isolated-tripod configuration 
as AAf(O) — (1/3) — 2^2- The spin chains in the centers 
of the tripod hexagons flip at the new critical magnetic 
field Be = 2J. While the magnetic fields increases fur- 
ther, the curve occurs at the 1/3-plateau and coincides 
with the curve obtained in the first approximation (see 
Fig. 12). 

The third approximation is to include configurations 
with isolated pairs of connected tripods. These configu- 
rations change the heights of the steps but give no new 
features in the magnetization curve. The probability of 
the pair of isolated tripod is P3 w 0.010. The entropy 
density increases up to 5*3 ~ 0.273. 

An isolated configuration of three tripods connected 
as a star is taken into account in the framework of the 




FIG. 3: An arrangement containing the tripod configuration. 
The tripod is shown in the solid line. The three white circles 
close to the center of the tripod have 4 spin-up (the black 
circles) and 2 spin-down (the white circles) nearest-neighbors. 




FIG. 4: Three tripods connected in the star arrangement. 
Three white circles close to the center of the star are sur- 
rounded by 5 spin-up and one spin-down nearest-neighbors. 



fourth approximation (see FigQl . The probability of this 
configuration is P4 k, 0.0035 and the entropy density in 
the fourth approximation is ^4 « 0.280. The key feature 
of this configuration is that there are chain spins that 
are surrounded by 5 spin-up and one spin-down nearest- 
neighbors. These chain spins remain stable up to the 
new critical magnetic field Be = 4J. To calculate the 
step of in the magnetization curve we have to calculate a 
number of various sequences of spin flips. The flnal value 
of the step is AAf(4) « 0.008 that is significantly smaller 
then AM(2) « 0.12. 

There is a variety of more complex configurations then 
the tripod. However, it should be mentioned that further 
approximations should change the heights of the steps 
but they can not cause new features in the magnetization 
curve. As it follows from Eq.Q there exist only four 
critical magnetic fields related to the four configurations 
of the nearest-neighbors shown in Fig|21 

The magnetization curve obtained in the fourth ap- 
proximation reproduces the key features of the experi- 
mental data at the very low sweep rate, namely, the four 
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equidistant steps in the magnetization curve. In contrast 
with our results, the third step in the experimental curve 
is much larger than the second one. This quantitative 
discrepancy can be eliminated in higher order approxi- 
mations that can be investigated both analytically, by 
the technique used in the present Letter, or numerically, 
applying Monte Carlo simulation . Approximations of 
higher orders are of importance for quantitative calcula- 
tions of the magnetization curve because the calculation 
convergence for the magnetic moment is slower than that 
for the configuration probability. 

Few questions on the magnetization of Ca3Co206 are 
still unclear. There were observed weak smeared features 
in the experimental curve at high magnetic fields above 
the last step. They can stem from the small fraction of 
Coll sites that are the high-spin state, because they in- 
crease the chain spin and cause higher critical fields. The 
sizeable hysteresis also draws attention in the experimen- 
tal curve. It depends drastically on the temperature and 
the magnetic field sweep rate. It should be mentioned 
that while the magnetic field sweeps down and crosses 
the highest critical field Be = 6J the chain spins flip 
down at random and the system occurs in a new state 
that is different from that at the sweeping-up process. 

In conclusion, we have developed a new model for 
the step-like magnetization of Ca3Co206 spin-chain com- 
pound. It can be applied also to other spin-chain 
compounds with the triangular lattice of chains, e.g. 
CaaCoRhOg. Due to the in-chain FM coupling between 
Co"^"*" ions at trigonal sites, C02O6 chains are considered 



at low temperatures as a large rigid spins with the strong 
Ising-like anisotropy. The AFM Ising model on the trian- 
gular lattice is applied to the system of rigid FM-ordered 
chains. The crucial point of the model is the supposition 
that the system is out of equilibrium, because the depen- 
dence of the magnetic moment on the magnetic field in 
the ground state of the AFM Ising model on the trian- 
gular lattice is smooth with the exception of the step at 
the zero magnetic field For the honeycomb AFM 

structure two steps were found in the theoretical mag- 
netization curve in excellent agreement with experimen- 
tal data at the intermediate temperatures. At higher 
approximations four equidistant steps were determined 
in accordance with experimental curves at the low tem- 
peratures and very low magnetic field sweep rate. The 
results obtained in the present Letter and the model of 
Ref.,2| can be regarded as two limiting cases. The first 
deals with the nonequilibrium interchain ordering assum- 
ing the chain spins to be rigid. An in-chain fragmentation 
and the quantum tunnelling of the magnetic moment of 
the fragments are investigated in the second case totally 
neglecting the interchain ordering. 
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